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anemometer positions are located along the bottom of the display. A data point above 
each anemometer position is deflected upwaxd by an amount proportional to tiie 
amplitude of the corresponding generator output voltage. A curve is passed through 
the vertically deflected points. The location of vortices with respect to individual 
anemometers can be deduced from the shape of the curve. GWVSS, however, is not 
operationally deployed because poles are not permitted near airport glide slopes for 
safety reasons. 

Another wake vortex detector of the prior art is a horizontally pointing Radar 
Acoustic Sensing System (KASS). Such a system has detected, tracke4 and measured 
the strength of wake vortices undo* all weather conditions. RASS radiates overlapping 
radar and acoustic beams. Because the air's index of refraction is a ftm.ction of 
density, the RASS acoustic beam, which consists of a spatial pattom of condensations 
and xarefrictions, generates in situ refractive index variations. Radar waves reflect 
from these index variations. As previously stated, radar reflections are strongest when 
the acoustic wavelength is half the radar wavelength. The reflections focus back onto 
the radar antenna because the radar and acoustic beams have overlapping spherical 
wavefronts. Radar reflections are Doppler shifled by an amount corresponding to the 
instantaneous speed of sound. Because a vortex's circular flow speeds up and slows 
down different parts of a RASS acoustic wave as it passes through the vortex, the 
radar Doppler spectrum is a mapping of the vortex's line-of-sight velocity distribution. 
Vortex circulation (vorticity) is deduced from the Doppler spectrum. RASS, however, 
emits a loud annoying sound, is costiy, and expensive to maintain. 

Wake vortices radiate audible sounds. These sounds have been recorded 6 and 
12 seconds after aircraft passage by a ground-based, vertically pointing 64-element 
microphone array. A system, having the acronym SOCRATES (Sensor for Optically 
Characterizing Ring-Eddy Atmospheric Turbulence Emanating Sound), for detecting 
this sound utilizes a laser-based, opto-acoustic technique to emulate a large array of 
microphones. SOCRATES has detected sound generated by wake vortices, but not 
ccmsistently. 
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Numerous tornado warning systems have been devised that rely on detecting 
sound generated by tornadoes. One such system, disclosed in U.S. Patent No. 
5,355,350 issued to H.E. Bass, et al, monitors outdoor noise between 180 and 420Hz. 
When the monitored sound intensity exceeds established thresholds, logic is applied 
to determine whether the level is increasing at a rate indicative of an ^proaching 
tornado and therefore whether an alarm should be sounded. This system, because 
sound in the selected passband is significantly attenuated by the atmosphere, detects 
tornadoes to a distance of only one-half mile and consequently, provides warning 
times of only 30 to 60 seconds. This method also does not provide range to the 
tornado. 

Another system for providing tornado warnings is disclosed in U.S. Patent No. 
6,097,296, issued to S. Garza, et al. In this system a narrow pass band acoustic filter 
is provided which is centered at one Hertz and designed to attenuate signals that are 
greater than and less than one hertz. In essence the filter has a near zero bandwidth. 
The wave emanating from this filter is clipped both at the peaks and valleys to obtain 
a wave that is substantially square. The square wave is coupled to a counter wherein 
the number of square wave cycles are counted. An alarm sounds when the counter 
reaches a predetermined number of cycles. 

This system is very difficult to implement. A filter that passes only one 
frequency, if possible, is extremely difficult to design. Further, a filter centered at IHz 
may not be optimum. The shape of a tormido's infrasound spectrum depends on 
internal wind speeds and tornado size, giving rise to infrasounds over a bandwidth that 
may not include one Hertz or may peak at some other frequency. Should the system 
be operable as described, it would be of limited utility since it does not provide range 
to the tornado and a concomitant warning time. 

CAT Detection : CAT is produced by high altitude atmospheric flows having 
speeds of 100 to 200 knots. CAT is often hundreds of miles long, tens of miles wide, 
and thousands of feet high and cannot be seen or avoided by pilots without prior 
knowledge of its location. 
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SUMMARY OF THE INVENTION 
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Figure 2A is a representation of the positioning of microphones at the middle 
marker of a glide slope. 

Figure 3 is a block diagram of a preferred embodiment of the invention for 
ground-based detection and location of wake vortices at an airport. 

Figure 3A is a dia^am indicating the position of microphone raws between 
parallel runways. 

Figure 3B is a diagram indicating the position of microphones in a column 
along the center line of a runway. 

Figure 4 is a block diagram of a preferred embodiment of the invention for 
ground-based detection and location of microbursts in the vicinity of airports. 

Figure 4A is a diagram indicating the positioning of parallel rows of 
microphones centered on tiie foot print of an arrival ghde slope. 

Figure 5 A and 5B are block diagrams of preferred embodiments of the 
invention for ground-based detection of tornadoes and for prediction of the time for 
a tornado to arrive at the sensing site. 

Figures 6 and 7 are block diagrams of a preferred mbodiment of the invention 
for early detection of CAT from an aircraft at high altitudes, the determination of CAT 
direction and distance from the aircraft, and the detection of microbursts from an 
aircraft on final approach. 

Figure 8 is a representation of relative positions of CAT and an aircraft:, and is 
useftil for the determination of range from the aircraft to the CAT region and the 
direction from the aircraft to the CAT region. 

DESCRIPTION OF THE PREFERRED EMBODIMENTS 

A flow is laminar or turbulent depending on the value of its Reynolds number 
Re, which is given by Re =UL/v, where X is a characteristic length of the flow, Uisa 
characteristic speed of the flow, and v is the fluid kinematic coefficient of viscosity. 
When the Reynolds number is below a critical value, viscosity smooths out lateral 
velocity variations, and laminar flow is favored. 
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made on top of Mount Washington. The 20 mph wind spectrum was obtained from 
sound recorded at Jones Beach (NY) in 2002. The 13 mph wind spectrum was 
obtained from recordings at JFK Airport. 

Figure 1 shows that the highest sound intensity produced by atmospheric wind 
occurs below 20Hz. If a spectrum's corner freqttency is defined by the intersection 
of the average fallofif slope of frequencies between 50Hz and 1000 Hz with a 
horizontal line approximating spectral an^litude between 4 and lOHz, it is evident 
from Figure 2 that the comer frequency shifts towards lower acoustic frequencies with 
increasing wind speed. This suggests that the population of the largest eddies in wind 
flow increases with increasing wind speed, resulting in a spectral shift of peak sound 
intensity to lower acoustic frequencies. The figure shows that each wind sound 
spectrum decreases about 6dB/octave in the audible band and that the intensity of 
sound below 20Hz is roughly proportional to the eightb power of the wind speed. 

The relatively high infrasound intensity at low frequencies for modest wind 
speeds shown in Figure 1 implies that an infrasound sensor can detect high speed 
hazardous turbulence at long ranges, for example, thunderstorms and hurricanes. 
However, both of these hazardous turbulences are heavily laden with moisture, 
making them easily detectable by radar. Infirasound sensing is particularly usefiil for 
long range detection of clear air hazardous turbulence, which is normally difficult to 
detect by other means. CAT, microbursts, tornadoes, and aircraft wake vortices are 
examples of clear air hazardous high speed flows that are difficult to detect by other 
means. 

CAT flows are free, while microburst, tornado, and wake vortex flows are 
generally free except near the ground. As mentioned above, turbulent flows within a 
few hundred feet of the ground generate slightly stronger sound than flow above the 
ground, enhancing detectability.. 

Wake Vortices: Sound generated by wake vortices were recorded by the 
inventor at JFK International Airport at a location 400 feet from the glide slope of a 
runway near the ghde slope middle marker. Sample sound spectra of tihese recordings 
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aircraft as obtained fi-om the laatched filter 3-4. This averaging is performed over a 
10 second time period that begins immediately after the gate 3-5 is qpened. The 
threshold is set by the threshold computer 3-7 at a value that is 20 dB lower than the 
measured average magnitude. This set threshold is coupled to the vortex detector 3-6, 
wherein signal magnitudes coupled firom the matched filter 3-4 are continuously 
compared to the set threshold. 

After the display unit 3-8 is activated, signal magnitudes coupled to the vortex 
detector 3-6 that exceed the set threshold are coupled thereto. The display on the 
display unit 3-8 has a graph-like appearance. Microphone positions are located along 
the bottom of the display and a vertical bar appears above each microphone position, 
the height of which is proportional to the corresponding measured infrasound 
magnitude. The approximate location of each vortex may be inferred by interpolating 
between vertical bar heights. When all five bar heights remain below a predetermined 
value for ^ pre- selected time period, hazardous vortices are deemed absent within 400 
feet of the runway center line around the middle marker. Display unit 3-8 is active 
only during the period of time that the gate 3-5 is open, it is turned off when the gate 
3-5closes. It should be recognized that the microphones need not be centered on the 
glide slope and their positioning is not limited to the middle marker. 

The FAA, in addition to imposing constraints between arrival times on a 
runway, also imposes time spacing constraints between aircraft arriving on parallel 
ruuM^ys less than 2500 feet part in order to prevent aircraft encounters with vortices 
that travel fi-om one runway to the other. Refer now Figure 3A. To monitor the 
existence of wake vortices between runways, parallel rows of microphones 3A-lmay 
be positioned between the two parallel runways 3A-2,3A-3, with the longitudinal 
center line 3A-4 of each row perpendicular to the runway direction 3A-5, that detect 
and track vortices crossing firom one runway to the other. Each row may have 3 
microphones. This number is not limiting, any number of microphones may be 
included in a row. Signal processing of the infrasound signals sensed by the 
microphones between the nmwajrs is similar to that previously described. 
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Opened. Opening the gate 3-5 activates the matched filter 3-4, the threshold computer 
3-7, and the display 3-8. If another aircraft is detected before the 3 minutes have 
passed, gate 3-5 stays open until fliree minutes have elapsed from the time the second 
aircraft was detected. 

Vortex detector 3-6 and the threshold computer 3-7 are coupled to receive the 
output signals from the matched filter 3-4. Vortex detector 3-6 continuously compares 
the magnitude of each signal to an ambient wind threshold which may be computed 
by threshold computer 3-7 from the infrasound generated by the ambient wind. The 
computation of the ambient wind threshold may be accomplished by averaging Ihe 
ambient wind infrasound magnitudes coupled from the matched filter 3-4 for each of 
the two microphones 3B-8,3B-9 fiirthermost from the departure threshold 3B- 
7during, for example, a 10 second time period that starts each time the gate 3-5 is 
triggered open and during which wake vortices have not yet been generated. It should 
be recognized, however, that a suitable fixed threshold for the vortex detector 3-6 may 
also be employed. 

When display 3-8 is activated by gate 3-5, signals from the vortex detector 3-6 
having magnitudes that exceed the second threshold are coupled to the display 3-8. 
As stated previously, the screen of the display 3-8 has a graph-like appearance. 
Microphone positions along the runway are located along the bottom of the display. 
A vertical bar appears above each microphone position, the height of which is 
proportional to the corresponding measured infrasound magnitude. During takeoff, 
the display shows the location of a departing aircraft as it accelerates down the 
runway. This is accomplished automatically since aircraft generated infirasound is 
coupled by successive microphones through the low pass filter 3-2, the matched filter 
3-4, and the vortex detector 3-6 to the display unit 3-8 as the aircraft makes it run 
down the runway until liftoff. Following liftoff, bar heights reveal the presence of 
wake vortices along the runway. When all bar heights remain below a predetermined 
level for a predetermined time, the departure corridor is deemed free of hazardous 
vortices. 
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When aircraft noise falls below the threshold set in engine noise detector 4-2, gate 4-4 
opens and signals firom the infrasound low pass filter 4-S are coupled to a coincidence 
detector 4-6, wherein each representative microphone signal magnitude is compared 
to a piedetertnined threshold having a value corresponding to the sound magnitude 
generated by, for exanq)le, a 44 knot ambient wind. 

When signals from all microphone representative signals exceed the 
predetermined threshold for predetermined time, which may be 30 seconds or more, 
a signal is coupled from the coincidence detector 4-6 to a microburst alert display unit 
4-7, whereon the location of all microphones receiving infrasound lhat exceed the 
threshold set in the coincidence detector 4-6 are displayed on a map of the airport 
glide slopes. 

Tornadoes '. A tornado some distance away is detected when its radiated 
infrasound arriving at a microphone exceeds a predetermined threshold. Lifrasound 
sensitive microphones similar to those described earlier may be used. Becai:^e cost 
is important, however, for this application, low cost electret microphones similar to 
Panasonic WM-034B may be used. Tests of these microphones show high pass 3dB 
points between 0.5 and O.S Hz. 

The SPL of infrasoimd emitted by a tornado having 120 knot internal winds is 
conservatively estimated to be 138dB. For a tornado with a diameter of 0.2 miles 
close to the ground, spreading loss reduces this value to an infrasound SPL of 98dB 
ten miles away, (138 -101og[10/,l]^)dB. A tornado sensing threshold based on this 
value, for example, corresponds to a 24 minute warning of a tornado travelling 25 
miles per hour towards the microphone. Based on the previous description of sotmd 
generated by turbulent flows, a tornado's spectrmn is expected to peak below 5Hz. 
For the detection of tornadoes, two embodiments will be described: The JEirst has 
simpler signal processing, while the second has a lower false alarm rate. Both 
embodiments provide tornado warning in excess of 20 minutes and also predict the 
time for a tornado to arrive at the microphone location. 
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predicted remaining time for the tornado to arrive at the microphone site calculated at 
time t{+i. The predicted time is; a positive number w^en r,.; - / is greater than 1. If a 
tornado bypasses the microphone site, the value of - 1 becomes less than 1 and the 
predicted time is a negative number, in. which case display 5-5 indicates that the 
tornado has passed. 

Refer now to Fig. 5B. Signals generated by the microphone are coupled to an 
analog low pass filter and analog to digital converter 5-9 wherein frequencies above 
lOOHz, for example, are converted to digital format at a 400Hz rate, for example. The 
digital signals at the output of the low pass filter and analog to digital converter5-9 are 
coupled to a fast Fourier transform (FFT) unit 5-10 wherein a 4096 point FFT is 
performed that is repeated every second. Based on a 500Hz A/D sampling rate, the 
waveform processed by the FFT algorithm each second is 8 seconds long. Only the 
outputs of frequency bios up through 15 Hz are retained. The FFT frequency 
resolution is about 0.12 Hz. The FFT 5-10 output is coupled to a calculator 5-11 
wherein the mean frequency, the 6dB bandwidfli, and the rms amplitude of the FFT 
bins within the 6dB bandwidth are calculated. These three computed quantities are 
coupled to a comparator 5-12 which determines whether the following three 
conditions are simultaneously satisfied over a period of 30 seconds: (1) the calculated 
rms bin amplitude within the 6dB bandwidth is greater than a predetermined value, (2) 
the calculated mean frequency is less than a predetermined value, and (3) the 
calculated 6dB bandwidth is greater than a predetermined value. When the preceding 
conditions are met, the comparator 5-12 provides a signal to open a gate 5-16 v/bich 
activates the tornado alarm and display 5-13. 

The data stream of calculated rms bin amplitudes is also coupled from the 
calculator 5-1 1 to a time delay unit 5-14 and to remaining time computer 5-15. The 
time delay unit 5-14 delays the data stream by a time AT which may be in the order of 
15 seconds and couples the time delayed data stream to the remaining time computer 
5-15. Usiag the delayed and undelayed data sfreams, the remaining time computer 5- 
15 calculates the *time remaining* before the tornado arrives at the microphone site 
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by 6dB to 93 dB to account for enhancement of aerodynamic wind sound g^erated 
near the ground (ground interaction). For a 37 mph surface wind speed, p== 1.293xl0- 
^ g/c^l^ V = 16.3x10^ cm/s, c= 3.3x10"* cm/s, and Af= v/c= 0.049. It follows from the 
free flow proportionality relation that CAT infrasound generated by an assumed 120 
knot wind flow at 30,000 feet is conservatively 135.9dB SPL, where p = 4.593x10-* 
g/cm^ V = 60x10^ cm/s, c= 3x10* cm/s, M=vA; = 0.2. For a 2000 foot high CAT 
flow, CAT infrasound is estimated to extend down to below IHz. 

The SPL power Py of sound arriving at an aircraft a distance y away from the 
center of a CAT flow of width A that radiates an SPL of P^, is &vea by die following 
formula: 

Py = {Pl + 10 log \1-Mcos0 +10 log iD/\y\)^] dB 

where: 0 is the angle of the observation point relative to the direction of the CAT 
flow. 

From the above, P^ is 135.9dB SPL. For a Macb 0.2 CAT flow (120 knots), the 
amplitude of the second term varies from +4.8dB to -4dB as a ftmction of 0. A value 
of -4dB is conservatively assumed. The third term quantifies sound spreading loss. 
For an aircraft 40 miles away from the center of a CAT flow 10 miles wide, the third 
term is 101og(10/40)^ = -12dB in which case Py = (135.9 - 4 -12) = 1 19.9dB. 

The above CAT SPL generates a microphone signal that is significantly greater 
than microphone self noise. It will now shown that CAT infrasound arrivmg at the 
microphones is much greater than that generated by the aircraft. Table 1 shows 
spectra over octave bands of aircraft noise radiated by a military jet aircraft at 1600 
feet altitude measured by the University of Kansas. These measurements are lOdB 
higher than comparable measurements made by the inventor on commercial jet aircraft 
arriving at JFK airport. 
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microphones are positioned approximately 50 feet in fiont of the eagines the SPL in 
each of the lower octave bands is increased by 101og(1600/150)2 = 20.6 dB to accoimt 
for the shorter distance traveled by the sound to the microphones and then subtracting 
an attenuation factor of 15.9 dB, resulting in a net increase of 4.7 dB in each octave 
band. With these assumptions, the total SPL of jet flow noise at the microphones in 
the combined 3-6 Hz, 6-12Hz and the 12-24Hz bands is (97 + 4.7 + 4.8) dB or 
106.5dB, which is conservatively 13.4 dB lower than the CAT SPL. 

For an aircraft travelling in a 40 knot ambient wind at 30,000 feet, the SPL of 
sound radiated by the ambient wind is (81.7 - 6)dB. The 6dB reduction accounts for 
the absence of wind sound enhancement produced by wind flow adjacent to the 
ground. The ambient wind SPL sensed by the onboard microphones varies from this 
value by lOlog |l-M;os^ 1*^ where ^ is the angle of aircraft travel relative to the 
direction of the wind. For the above example, the SPL varies with wind direction from 
+1.5 dB to -1.4dB. This variation is smoothed out by the 360 degree sound 
acceptance angle of each microphone. Combining all of the above, a 40 knot ambient 
wind generates an SPL about 39dB lower than the CAT SPL. 

The Doppler effect produced by forward aircraft motion transforms each 
arrivkig sound frequency /into Ff , w^ere F = (IWcos^ ), M is aircraft Mach 
number, and d is the direction of sound arrival relative to the direction of aircraft 
travel. The amplitude of sound at each Doppler transformed frequency is multiplied 
by the same factor. Since ambient wind sound arrives at the microphones from all 
directions, the frequency and amplitude Doppler effects on the sensed ambient wind 
sound average out. 

For the case in which CAT sound arrives in the aircraft's forward hemisphere, 
the value of F varies between 1 and 1.8 for a Mach 0.8 aircraft. For F =1.8, a CAT 
sound spectral peak around 0.5Hz, for example, is transformed into 0.9Hz with a 
2.6dB SPL increase. Thus, the Doppler effect introduces a relatively small spectral 
shift and increases CAT infrasound detectability. 
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Refer now to Figure 6. Signals from microphones positioned on an aircraft 6- 
lA and 6-lB are coiqiled, respectively, to low pass filters and analog to digital 
converters 6-2A and 6-2B where ttiey are (1) analog low pass filtered to eliminate 
signals above lOOHz, for example, (2) sampled at a 2048Hz rate, for example, which 
satisfies the Nyquist sampling criterion, and (3) converted to digital format The 
digitized signals are coupled from the filters and converters 6-2A and 6-2B, 
respectively, to FFT imits 6-3 A and 6-3B, each of which may perform a 2048 point 
Fast Fourier Transform (FFT) every second. FFT output bin frequencies may be 
spaced IHz apart Only the outputs from FFT bins having frequencies in the 1 to 25 
Hz band are retained, though other bands may be utilized. 

FFT bin amplitudes and FFT bin phases affor frequencies in the l-25Hz band 
are coupled from the FFT unit 6-3A to a comparator 6-4A, along with their frequency 
identification. Similarly, FFT bin amplitudes and phases in the l-25Hz band are 
coupled from FFT unit 6-3B to comparator 6-4B, along witii their frequency 
identification. Comparator 6-4A compares each bin amplitude a^to a predetermined 
threshold. When a bin amplitude is less than the threshold, the bin's amplitude and 
the corresponding bin phase are eliminated from the retained set of FFT bin outputs. 
This eliminates signal outputs from bins having low S/N. Comparator 6-4B performs 
a similar operation. 

Conq)aratar 6-4A couples retained bin amphtodes with their bin frequency 
identification to a nns amplitude and average frequency computer 6-5A and to a 
weighted angle ofiFboresight computer 6-9, and also couples corresponding FFT bin 
phases with their bin frequency identification to a phase difference computer 6-7. 
Similarly, comparator 6-4B couples retained bin amplitudes to a rms amplitude and 
average frequency computer 6-5B and to the weighted angle off boresight computer 
6-9, and corresponding FFT bin phases to the phase difference computer 6-7. Rms 
amplitude and average frequency computer 6-5A calculates (1) the rms sum of 
retained bin anq)litudes, denoted by A, and (2) the bin ampUtude weighted mean 
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fiequency, denoted by F^. Rms amplitude and average ftequency computer 6-5B 
perfonns likewise, generating ims sum B and mean frequency Fg. Amplitude A and 
frequency Fa are coupled to a second comparator 6-6A. Amplitude^ is also coupled 
to rms averager 7-1 in Fig.7. Similarly, amplitude B and frequency Fg are coupled to 
comparator 6-6B. Amplitude B is also coupled to llie rms averager 7-1 in Fig.7. 

When signals coupled from FFT units 6-3 A and 6-3B to con^arators 6-4A,6-4B 
exceed the predetermined thresholds, the frequencies F^ and F^ established, 
respectively, by the computers 6-5A and 6-5B are measures of the largest scale eddies 
in the turbulence generating the detected sound. Should hazardous CAT not be 
present, the rms sum and the average frequency determined in each microphone 
channel are representative of the ambient wind. For this case, bin amplitudes provided 
by the FFT units will be lower than the predetermined thresholds and will not be 
coupled to the rms ampUtude and average frequency computers. As an aircraft 
approaches a CAT region, the FFT bin amplitudes will shift from (hat characterizing 
ambient wind towards ampUtudes that characterize CAT and the processing described 
above will continue. The FFT bin amplitudes will effectively reach the CAT 
amplitudes i^en the received CAT SPL is about 6dB greater than the ambient wind 
SPL. 

As (he aircraft approaches closer to the CAT region, the rms sum m each 
channel will continue to increase. Comparator s 6-6Acompares the measured mean 
frequencies coupled from the computer 6-5 A to a predetermined threshold frequency, 
which may be set at ^proximately 5 Hz. If the mean frequency falls below the 
threshold frequency, comparator 6-6A then compares the rms sum to a predetermined 
amplitude threshold. If the sum exceeds the threshold, a CAT alarm is sounded in 
channel A. 6-6B performs similarly. When a CAT alarm is set off in both channels, 
gate 6-10 opens. 

As mentioned above, comparators 6-4A and 6-4B each co\q)le retained bin 
phase values to the phase difference determinator 6-7 with their corresponding bin 
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frequency identification. Phase difference detenninator 6-7 computes the electrical 
phase difference 0f between each pair of retained phase values having the same 
frequency identification, when both phase values are present. Electrical phase 
differences <2^and associated bin frequency identifications are coupled to a converts 
6-8 vsdlierein tiie phase difference is converted to the angle off boresight of the CAT 
sound arrival. For an assumed microphone separation of 10 feet, the maximum value 
of ^for any frequency up to 25Hz is less than 30 electrical degrees (at 30,000 feet). 
Each electrical phase difference <Z^is converted to an estimate of CAT sound arrival 
angle Of (measured from aircraft boresight) using the following formula: ^ = cos-*(^ 
c„ I14D\ where c„ is the speed of sound at aircraft altitude,/is frequency, and £> is the 
horizontal separation between the two microphones. 

Arrival angle estimate Of axe coiqjled &om the converter 6-8 to a weighted 
angle off boresight computer 6-9 to which each corresponding identifying frequency 
is also coupled. Computer 6-9 also receives from each of comparator 6-4A and 6-4B 
retained bin signal amplitudes and Ofg, respectively. Computer 6-9 calculates the 
average of 0,51 and Oyg, denoted by Afi for each retained frequency and then calculates 
a signal strength weighted estimate of the arrival angle, denoted by 0 , using the 
following formula: 0 = £Af Of I EAf. Because the microphones have no directivity, 
the computed value of 0 is ambiguous. Depending on whether the sound was received 
in the forward or back hemisphere, the correct value is either 0 or [180- ^, 
respectively. This ambiguity is resolved by the processing indicated in Figure 7, with 
reference to the geometry shown in Figure 8. 

Referring Figures 7 and 8, an aircraft is shown at position 1 when a CAT alarm 
is announced in botii microphone channels, and at position 2 at a time <4r seconds 
later. The aircraft travels a distance X between the two positions at an aircraft ground 
speed of Kknots, so that X-VAT. Position 1 is a distance from the CAT region 
at which time the rms sum signal of the bin signal amplitudes at the computer 6-5 A 
is Ai and the rms sum signal of the bin signal ampUtudes at the computer 6-5B is Bi. 
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couples this value to the distance to CAT computer 7-5 and to the gate 6-10. The 
distance to CAT computer 7-5 then calculates using the value [180 -6J and couples 
this result to the gate 6-10. As a result the range and angle to CAT are coupled 
through gate 6-10 to provides updated distances and unambiguous boresight angles 
to the CAT region while gate 6-10 remains open. 

When the aircraft is traveling almost parallel to the region, 1 6 1 will be close to 
90** in which case the distance to CAT region will change very slowly, giving 
inaccurate results for computed distance to the CAT region. Hence, when |r-l| is 
equal to or less than 0.05, the distance prediction previously displayed is not updated 
until |r-l I exceeds 0.05. 

The embodiment for detecting CAT and measuring distance and direction to a 
CAT region may also be utilized onboard an aircraft to detect microbursts during final 
approach and to measure the distance and direction to the microburst region. Since 
microburst wind speeds are lower than CAT flow speeds, the radiated aerodynamic 
sound intensity will be smaller, resulting in shorter microburst detection ranges. 
Aircraft landing speeds, however, are approximately one third of the speeds at high 
altitude, so that microburst warning times may be only slightly less than CAT warning 
times. 

It should be recognized that frequencies, times, and other parameters mentioned 
in the description of the preferred embodiments are given by way of example, whether 
or not so specified, and that other such values may be used to provide the desired 
results. 

While the invention has been described in its preferred embodiments, it is to be 
understood that the words which have been used are words of description rather than 
of limitation and frequencies, times, and other parameters are given by way of 
example, whether or not so specified. Thus the scope of the invention is determined 
by the appended claims and that changes may be made within the purview thereof 
without departure from the true scope and spirit of the invention in its broader aspects. 
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I claim: 

1 A method for detecting atmospheric disturbances including the steps of: 
providing infirasoimd frequency magnitudes of received noise spectra; 
comparing said infrasound frequency magnitudes to a computed infrasound 

threshold; and 

determining locations of said atmospheric disturbances with the utilization of 
infrasound frequency magnitudes that exceed said threshold. 

2. A method for detecting atmospheric disturbances in accordance with claim 1 
wherein said providing step includes the steps of; 

extracting noise at frequencies below a specified frequency from said received 
noise spectra to provide an extracted noise spectra; 

fatering said extracted noise spectra to obtain infrasound at frequencies below 
a predetermined infrasound frequency; and 

detecting magnitudes of infrasound frequencies below said predetermined 

infrasound frequency. 

3. A method for detecting atmospheric disturbances in accordance with claim 2 
wherein said extracting step includes the step of activating said filtering step when 
magnitudes of said extracted noise spectra exceed a preselected threshold. 
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